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Abstract We study Lagrangian statistics of the magnitudes of velocity and pres- 
sure gradients in isotropic turbulence by quantifying their correlation functions and 
their characteristic time scales. In a recent work [1], it has been found that the La- 
grangian time-correlations of the velocity and pressure gradient tensor and vector el- 
ements scale with the locally-defined Kolmogorov time scale, defined from the box- 
averaged dissipation-rate (er) and viscosity (u), according to ^ — \/v l^r- In this 
work, we study the Lagrangian time-correlations of the absolute values of velocity 
and pressure gradients. It has long been known that such correlations display longer 
memories into the inertial-range as well as possible intermittency effects. We explore 
the appropriate temporal scales with the aim to achieve collapse of the correlation 
functions. The data used in this study are sampled from the web-services accessible 
public turbulence database (http : //turbulence .pha. jhu. edu). The database archives 
a 1024^ (space-|-time) pseudo-spectral direct numerical simulation of forced isotropic 
turbulence with Taylor-scale Reynolds number Re\ = 433, and supports spatial differ- 
entiation and spatial/temporal interpolation inside the database. The analysis shows 
that the temporal evolution of the auto-correlations of the absolute values are deter- 
mined not by the local Kolmogorov time-scale but by the local eddy-turnover time scale 

2 /o —1/3 

defined as re,r = r ' . However, considerable scatter remains and appears to be 

reduced only after a further (intermittency) correction factor of the form of (r/L)^ is 
introduced where L is the turbulence integral scale. The exponent x varies for different 
variables. The collapse of the correlation functions for absolute values is, however, less 
satisfactory than the collapse observed for the more rapidly decaying strain-rate tensor 
element correlation functions. 
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1 Introduction 

The study of turbulence from a Lagrangian viewpoint has a long history, with the 
earliest works of Taylor [2] and Richardson |3 both pre-dating Kolmogorov [4]. The 
Kolmogorov 1941 theory used the constancy of the globally-averaged dissipation-rate 
(e) across scales to deduce, among others, the scaling properties of the wavenumber 
spectrum of kinetic energy. The 1941 theory was extended to account for intermittency 
by the introduction of the so-called refined Kolmogorov similarity hypothesis (RKSH) 
[5] in 1962. In this K62 extension of the theory, conditional statistics, based on the 
dissipation rate averaged in some particular subregion of the fiow, acquires a central 
role. The local dissipation rate, usually denoted by e^, is defined according to 



where V is the volume of the subregion 7?.r(x) (e.g. a box or a sphere) of size r centered 
at X, V is the kinematic viscosity of the fluid, and Sij is the strain-rate tensor. One 
of the main predictions of the RKSH relates to the longitudinal velocity increment 
at scale r, defined as 5rU — [ui(x -I- r) — Ui{x)]{ri/r). The RKSH states that in the 
inertial range of turbulence the statistics of 5rU depend on r and e^- Therefore, from 
dimensional analysis, various moments of SrU conditioned upon a fixed value of er will 
scale as (SrU^ler) ~ Cp{rerY^^ , essentially following Kolmogorov's 1941 postulate, but 
locally. Anomalous scaling results from the additional global averaging and anomalous 
scaling behavior of moments of er. The existing literature to validate RKSH has focused 
mainly on velocity increments [6][7ll8ll9l ll0llll) or acceleration [12) in which the analysis 
can be performed on single snapshot measurements of the turbulent flow, i.e. based on 
a relatively 'static' point of view of the flow. 

In order to examine the RKSH in more depth, one would also like to inquire about 
its dynamical origin, specifically its role in the time evolution of the local structure of 
turbulence. It is well recognized that the dynamics of turbulence is best understood in a 
Lagrangian frame of reference, i.e. following fiuid particles. Recent years have witnessed 
a strong revival of interest in Lagrangian statistics in turbulence. For reviews, see |13l 
I14j . The dynamics of turbulence following fluid particles also plays a central role in the 
PDF modeling framework developed over the past two decades by S. Pope, starting 
with his seminal 1985 paper |15j . 

Besides which 'frame' to use in the description of the dynamics, it is also important 
to select variables of interest that convey rich information about the flow. In recent 
years, there has been growing attention placed in the dynamical evolution of the ve- 
locity gradient tensor A (Ajj = dui/dxj). This is due to the fact that A provides rich 
information about the topological and statistical properties of small-scale structure in 
turbulence. Pioneering studies of the Lagrangian structure and stochastic modeling of 
A are described in Refs. [16] and [17], respectively. The Lagrangian time evolution of 
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A can be obtained by taking gradient of the NS equation |18j : 

^ = -^^H-r.^^. (2) 

As usual, d/dt stands for Lagrangian material derivative, p is the pressure divided 
by the density of the fluid, and the second and third terms on the right-hand-side 
of this equation are the pressure Hessian tensor and viscous term respectively. Ne- 
glecting viscous effects and assuming the pressure Hessian isotropic lead to a closed 
formulation known as the Restricted-Euler (RE) dynamics [18], [19]. With analytically 
treatable solutions for the full tensor-level, the RE system provides a fruitful starting 
point for small structure modeling although there exist serious deficiencies in the RE 
dynamics, especially since it predicts nonphysical finite-time singularities [TU]. Models 
have been developed to mimic the regularization features of the neglected pressure 
Hessian and viscous terms. Efforts include a stochastic model in which the nonlinear 
term is modified to yield log-normal statistics of the dissipation [T7], a linear damping 
model for viscous term [2D], a tetrad model [5T] for pressure Hessian closure, a viscous 
diffusion closure [22], a new stochastic dynamic model, so-called Recent Fluid Defor- 
mation closure, for both viscous and pressure Hessian terms [23, and a multi-scale 
model which includes energy exchange between scales [2J . The study of the temporal 
auto-correlation structure of various quantities associated with A assists in the further 
developments and improvements of such models. 

It has recently been confirmed PP that when averaging over the entire domain, 
auto-correlation functions of velocity gradient tensor elements decay on timescales on 
the order of the mean Kolmogorov turnover time scale. This time scale is computed 
from the globally averaged rate of dissipation and viscosity. However, when performing 
the analysis in different subregions of the fiow, turbulence intermittency was found 
to lead to large spatial variability in the decay time scales. Remarkably, excellent 
collapse of the auto-correlation functions is recovered when using the 'local Kolmogorov 
time-scale' defined using the locally, rather than the globally, averaged dissipation-rate 
(r^ = ^Jvjer). This is an additional new evidence for the validity of Kolmogorov's 
Refined Similarity Hypothesis, but from a Lagrangian viewpoint that provides a natural 
frame to describe the dynamical time evolution of turbulence. 

In this paper, we study Lagrangian time-correlations of scalar measures (such as 
magnitudes) of velocity and pressure gradients and explore whether there is further 
evidence for the Lagrangian RKSH for these variables. Lagrangian correlation functions 
of square strain- and rotation-rate have aheady been studied in prior work [25ll26j . In 
examining the scaling of the magnitudes of the velocity gradient tensor and pressure 
gradient, the behavior of correlation functions will be shown here to be much more 
complex than for the tensor or vector elements themselves. Nontrivial dependencies 
on local-length scale r and local dissipation-rate (er) are observed, and these require 
more detailed study. The present paper is devoted to such a study, based on analysis 
of Lagrangian data. 

Lagrangian data can be extracted from direct numerical simulation (DNS) of NS 
equations with relative ease. The first such effort traces back to Riley and Patterson 
[27j . The rapid development in computing power over the past few decades has spurred 
vast amount of such numerical investigations at increasing Reynolds numbers. For 
relevant reviews, see [28lll3|[T4ll29j and references therein. A new way to exploit large 
databases in turbulence has been recently proposed [30]. This approach is based on 
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web-services that allow public access to turbulence DNS databases that store not only 
snapshots of 3D distributions but also the entire pre-computed time history. Using this 
public turbulence database, here we study the Lagrangian time evolution of velocity 
and pressure gradients and their magnitudes in isotropic turbulence. 

The remainder of this paper is organized as follows. Section [2] describes the public 
turbulence database and the numerical approach we use to perform the Lagrangian 
analysis. Results on the time evolution of auto-correlations of velocity and pressure 
gradient magnitudes are presented in Section |3] We conclude in Section |4] with a short 
discussion. 

2 JHU public turbulence database and particle-tracking approach 

The DNS data of a forced isotropic turbulence archived in the JHU public database 
system are from a pseudo-spectral parallel computation of the forced NS equations 
in a [0, 27r]'^ domain, at a Taylor-microscale Reynolds number of Re^ — 433 [30| . 
The database contains of output on 1024^^ spatial points and 1024 time samples (ev- 
ery tenth DNS time-step is stored) spanning about one large-scale eddy turnover time. 
The domain- wide averaged dissipation-rate (< e >) and the corresponding Kolmogorov 
time scale (tx) are 0.092 and 0.045 respectively, in the units of the simulation. The 
turbulence integral scale is L = 1.376. Some data processing functionalities such as 
spatial differentiation, and spatial and temporal interpolations are provided directly 
inside the database. This feature not only reduces data download cost but also allows 
users to obtain desired quantities at arbitrary locations and times. The whole database 
results in a 27 Terabyte storage size. The 1024^ space-|-time history of turbulence is 
publicly accessible through a web-service interface which serves as a bridge to connect 
user requests with the database nodes. Users may write and execute analysis programs 
using prevailing languages C, Fortran, or Matlab on their host computers such as desk- 
tops or laptops, while the programs request desired outputs from the database through 
GetFunctions (subroutine-like calls) over the Internet. Currently, eight GetFunctions 
listed in Table [l] for velocity and pressure along with their derivatives and force are 
available. With these call functions, users can retrieve quantities simultaneously for 
large amounts of locations and time (within the stored time frames) without expensive 
memory and time costs. The details of the DNS data and JHU turbulence database 
can be found in a previous publication [30]. The instructions and sample codes in C, 
Fortran, and Matlab are available at http://turbulence.pha.jhu.edu. 

We employ the particle-tracking algorithm of Ref. [31] to extract Lagrangian in- 
formation along many particle trajectories simultaneously. Each particle is tagged and 
randomly assigned an initial position. Let x'^(y,f) and u^{y,t) denote the position 
and velocity at time t of the fluid particle originating from position y at initial time to 
with the superscript -|- representing Lagrangian quantities following the fluid particle. 
Each particle is tracked by numerically integrating 

^^ = u+(y,t) (3) 

where the Lagrangian velocity u^(y,t) is replaced by the Eulerian velocity u(x, t) 
where the particle is located, namely u^{y,t) = u{x'^ {y , t) , t) . 

The particle displacement between two successive time instants tn and tfi-^i{^ — 
tn + St) is obtained through an integral of Eq. ((Sjl using a second-order Runge-Kutta 
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Function name 


Spatial difT. 


Spatial int. 


Temporal int. 


Outputs 


Get Velocity 




Noint, Lag4,6,8 


Noint, PCHIP 


Ui 


GetVelocityAndPrcssurc 




Noint, Lag4,6,8 


Noint, PCHIP 


Ui,P 


GetVelocityGradient 


FD4,6,8 


NoInt, Lag4,6,8 


Noint, PCHIP 


dxj 


GetPressureGradient 


FD4,6,8 


Noint, Lag4,6,8 


Noint, PCHIP 


(Jp 

dxi 


GetVelocityHessian 


FD4,6,8 


Noint, Lag4,6,8 


Noint, PCHIP 


dx^dxj 


GetPressurcHessian 


FD4,6,8 


Noint, Lag4,6,8 


Noint, PCHIP 


dx-i dxj 


GetVelocityLaplacian 


FD4,6,8 


Noint, Lag4,6,8 


Noint, PCHIP 


dxj dx j 


GetForce 




Noint, Lag4,6,8 


Noint, PCHIP 


h 



Table 1 Subroutine-like call functions, diff: difTeretiation; int: interpolation; Noint: no inter- 
polation; FD: Centered finite difference, options for 4th-, 6th-, and 8th-order accuracies; Lag: 
Lagrangian polynomial interpolation, options for 4th-, 6th-, and 8th-order accuracies; PCHIP: 
Piecewise cubic Hermite interpolation. 



method. At time tn for a particle located at x^(y,tn), the predictor step yields an 
estimate x* = x^(y,tn) + 5t u~''(y,t) for the destination position x^(y, t„-|_i). The 
corrector step then gives the particle position at tn+i: x^(y,t„+i) = x^(y,t„) + 
5t [u^(y,t„) + u~'"(x*, t„+i)]/2. It is proved that the time-stepping error is of order 
{Sty^ over one time step [31]. In general, accurate spatial and time interpolations are 
crucial to obtain the fluid velocities while tracking particles along their trajectories. 
In JHU turbulence database, these operations have been built in with optional orders 
of accuracy. We use flags FD4Lag4, LagS, and PHCIP (explained in the caption of 
Table [1} for the calls to specify spatial differentiation, spatial interpolation, and time 
interpolation. 



3 Lagrangian time correlations of gradient magnitudes 

As mentioned above, the velocity gradient tensor A has received considerable atten- 
tion in recent years. It provides a rich characterization of the topological and statistical 
properties of the small-scale structure in turbulence in the viscous range. The antisym- 
metric part of the velocity gradient tensor, i.e., fl^j = {A^j — Aji)/2, is the rate of 
rotation describing the vortex structure and dynamics. Whereas the symmetric part 
of A, the strain-rate tensor, defined as Sij = {A^j + Aji)/2, represents the strength 
and directions of fluid deformation rates. The dynamic evolution of A is given by Eq. 
([2]). Here we study Lagrangian auto-correlations for the absolute values of velocity- 
derivative tensor using the magnitudes of the strain-rate tensor S and rotation-rate 
tensor fl. The magnitudes are defined here using the square invariant according to 
|S| = ^ Sij Sij and |f2| = ^ Slij f2ij (we apply Einstein notation for repeating in- 
dexes unless indicated otherwise). We also study the magnitude of pressure gradient 
(approximately similar to the acceleration magnitude) , defined as | y p| = v^Vp~Vp. 
The Lagrangian time correlation of these scalar quantities is defined as usual: 

^ </(^o)/fa+-)) (4) 

^ V{mr}-{f{to+rr) 

where r is the time-lag along Lagrangian trajectories, / can be / = |S|, / — \fl\, or 
/ — 1 V p! s-s the case may be, and {• • •) may represent ensemble or global volume 
averaging for homogeneous turbulence. 
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In order to study effects of intermittency, which in turbulence is characterized by lo- 
cal regions displaying different levels of turbulence activity, we also compute conditional 
correlation functions based on fluid particles that originate from various subregions of 
the flow domain. The subregions are characterized by the local dissipation-rate er de- 
fined in Eq. ([T]) and the length-scale r. With this context, the global average in Eq.Q 
is replaced by the conditional average, i.e 

^ ^if{t^?\^r)-{f{to+rf\er) 

In the conditional average the initial position of particles contributing to the av- 
erage at time to are sampled from several local boxes of size r that have a prescribed 
locally averaged dissipation-rate tr- For practical reasons, a finite range of values of er 
must be considered, i.e. we use sampling in bins of tr values. For each bin there are sev- 
eral such local cubes for which their er falls in a prescribed range. Besides varying the 
bin location, we also consider four length scales, r = ZArix, GSrjj^, ISGrjj^, 272?7/f . They 
correspond to 16-, 32-, 64-, and 128 grid-point cubes, respectively. Each has four asso- 
ciated er bins. An additional set of 64-cube cases (r — 13677;^) with five bins has been 
studied. This set was also considered in [T] for the tensor element-based correlations. 
The specifications of scale and bin values are listed in Table [21 



cube size r 


bin index 


particle/cube 


cube/er 


Er range 


nominal 


34r?A' 


1 


50 


120 


0.0056 ^ 0.0093 


0.0076 




2 


50 


120 


0.017 0.020 


0.019 




3 


50 


120 


0.089 ^ 0.096 


0.092 




4 


50 


120 


0.14 ^ 0.16 


0.15 




1 


100 


60 


0.0074 ^ 0.011 


0.0096 




2 


100 


60 


0.0241 0.028 


0.019 




3 


100 


60 


0.089 ^ 0.096 


0.093 




4 


100 


60 


0.14 ^ 0.16 


0.15 


136r?A-I 


1 


500 


12 


0.017 0.020 


0.018 




2 


500 


12 


0.024 ^ 0.028 


0.026 




3 


500 


12 


0.089 ^ 0.096 


0.093 




4 


500 


12 


0.10 0.12 


0.11 


136r/x-n 


1 


4000 


2 


0.033 ^ 0.041 


0.037 




2 


4000 


3 


0.045 0.046 


0.046 




3 


4000 


2 


0.065 ^ 0.066 


0.066 




4 


4000 


2 


0.083 ^ 0.088 


0.086 




5 


4000 


3 


0.12 ^ 0.16 


0.14 




1 


1000 


6 


0.030 ^ 0.037 


0.034 




2 


1000 


6 


0.048 ^ 0.060 


0.052 




3 


1000 


6 


0.081 ^ 0.10 


0.094 




4 


1000 


6 


0.11 ^ 0.15 


0.12 



Table 2 Indication of different cases for calculation of conditional correlation functions. Shown 
for local subregion (cube) length size, bin index, number of particles in each cube, number of 
cubes in each bin, and range and nominal of local dissipation for each bin. 



In Fig.[T]we show PDFs of er for the four length scales considered (left plot) and 12 
representative 64-cubes placed inside the 1024^ domain, with 50 sample fluid particle 
trajectories emanating from each and progressing during a time equal to 27tk (right 
plot). The required averages are taken over all the trajectories as well as over several 
cubes for which er is in a bin's prescribed range. 
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Fig. 1 Left plot: PDFs of locally-averaged dissipation-rates with four different local length 
scales; Right plot: Sample particle trajectories starting from 12 randomly selected 64-cubes 
characterized by local dissipation-rate tr at the initial time corresponding to case ISGtjk — // 
in Table 12 



Lagrangian and conditional Lagrangian auto-correlation functions for tensors A, S, 
ft, or vector VP have been studied in our previous work pil32). These tensor and vector 
time correlation functions are computed through expressions like (C^^ (to)C.y (fo + r)) 
and {Cij{to)Cij{to + T)\er) or {Gi{to)Gi{to + r)) and {Gi(to)Gi(to + r)\er) on tensor 
or vector element level for global and conditional correlation functions, respectively. 
Here we present these measurements in Fig. [2] for three variables including all the cases 
listed in Table [2] (in [T] only results for S were shown for the same cases). The evolution 
time is scaled by tk and TK,r where 




The main observations can be summarized as follows. First, the rotation-rate displays 
significantly longer time memory than the strain-rate. After about Gr^^^j the strain- 
rate's correlation is essentially zero, whereas it is still near 0.5 for the rotation-rate. 
We found that this trend holds true even if coherent vortex structures are excluded 
from the analysis. Second, the temporal auto-correlation functions scatter significantly 
when the time lag is scaled by tx (left column) but collapse well when scaled by 

r (right column) . This behavior demonstrates that the dynamics of flow variables 
such as velocity and pressure gradients following fluid particles depends upon the local 
dissipation-rate (e^) rather than the global one ((e)) which, as argued in [T], provides 
new evidence for the validity of Kolmogorov's refined similarity hypothesis form a 
Lagrangian viewpoint. 

In what follows, we study conditional Lagrangian time correlations for the absolute 
values of these tensors and vector. 

Figure [3] shows the similar plots to Fig. [2] but for the the absolute values of tensors 
S (top row) and fl (middle row) and vector Vp (bottom row). It is quite clear that for 
all variables, especially for the strain-rate and pressure gradient, the correlations decay 
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Fig. 2 Lagrangian time correlations of strain- and rotation-rate tensors and pressure gradient 
vector for all the cases in Table [2] Open squares are for global average over randomly located 
particles in the whole domain, whereas different lines correspond to subregions of the flow 
characterized by different e^s. Time-lag is normalized using the global Kolmogorov time scale 
Tx (left column) and the local time-scale rx,r (right column). 



much more slowly for the magnitudes as compared to the tensor or vector elements. 
Similarly slow decay had been observed for the square of these variables in [25ll26j . 
Moreover, and unlike the tensor- or vector-based Lagrangian auto-correlations, poor 
collapse is seen when the time lag r is scaled by the local Kolmogorov time t^ j.. Such 
scaling appears to work only for the viscous time-scale range near the origin of the 
curves (r < tk ,.). For the inertial range, the curves scatter significantly even when 
scaled by the local r- 

Since there remain significant correlation scatters even after long time delays, we 
explore the use of other time-scales to express time. The characteristic time-scale that 
is believed to be relevant in the inertial range is the eddy-turnover scale appropriate 
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for eddies of size r. Its global and local values are defined according to 

r2/3/ \-l/3 2/3 -1/3 ,„N 

Te = L ' (e) ' , Te^r = r ' er '. (7) 

Fig. m shows the conditional Lagrangian time correlations of absolute values of S 
(top row), f2 (middle row), and VP (bottom row) with time normalized by re (left col- 
umn) and re,r (right column). These results are for a single length-scale corresponding 
to the case of r = 136?7x (case II) in Table [5] When the time lag r is scaled by Te, 
there are noticeable differences in the results depending on er - For the three variables, 
larger values of er (i.e. in regions of more intense turbulence activity corresponding to 
smaller local eddy turn-over time) are associated with faster correlation decay. When 
the time lag is scaled by the local time-scale Te,r, the curves collapse better than with 
the global value. This provides some evidence for a Lagrangian RKSH also at inertial- 
range 'eddy-turnover' scales. 
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The next question is whether good collapse also occurs for different length scales. 
Fig.[5]plots the conditional auto-correlations of and |Vp| with approximately 

the same er corresponding to bin No. 3 (see Table [21, but at different length scales. 
In the left column, the time lag is scaled by the local eddy time Te,r- It is seen clearly 
that the correlation functions decay differently at different length scales, which implies 
that the normalization of time with Te,r does not account for the differences. 

Intermittency in turbulence is often known to connect the inertial range dynamics 
with the ratio of length-scale to the integral scale, i.e. the level of intermittency is 
related to "how far" the scale is from its original starting point at the large scales 
during the cascade. Often such effects are parameterized by factors of the form [r/L)^, 
where x is an appropriate intermittency exponent for the correction. We determine the 
exponents x empirically (see below) to obtain improved collapse. The right column in 
Fig. [5] shows the correlation functions with time now normalized by the intermittency 
corrected time scale Te,r{r / L)^ . As can be observed, improved collapse can thus be 
obtained by using an intermittency correction. 

In order to determine the exponents empirically, we set a threshold on the correla- 
tion function value. A value Pfij^j-i) — 1/2 is used, which defines the time-scale T^j^. 
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Fig. 5 Time correlations of \S\, and |Vp| for the e,- value in bin No. 3 (see Table O 
vs. r/Tg^r (left column) and r/T^^r ■ (x)^ (right column, (b):x = —0.3 , (d):x = —0.49, and 
(f)-X = —0.55). All are for for the cases in bin No. 3 with e 0.092, which also close to value 
of global-averaged dissipation-rate. 



For each case, the value of ri^2 correspond to — 1/2. A log-log plot of Ti/2/'''e,r 
versus r/L should have slope Xi if a power-law intermittency correction is appropri- 
ate. In Fig. [6] such plots are presented for each variable |S|, |n|, |Vp|, from top to 
bottom. Exponents are obtained by fitting straight lines through the data as shown 
by the dashed lines in the plots. The corresponding exponents of x a^re X|s| = —0.3, 
X\Q\ = —0.49, and X|Vp| ~ —0.55. These are the values used to scale the results shown 
in right column of Fig. \S\ We have not yet succeeded in relating the values of x to the 
multifractal theory of turbulence. Note that these values are significantly larger than 
what is typically obtained from multifractal corrections. 



All the cases are plotted jointly in Fig. [7] With global time scaling of Te (left 
column) , the curves scatter significantly refiecting clear intermittency in the fiow. When 
both Te,r and (r/L)^ are considered in the scaling of time lag (right column), the curves 
collapse reasonably well. Some scatter remains, however. 
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Log,o(r/L) 



Fig. 6 Plots of characteristic decay time r-^ji (measured as the 1/2 point in the correlation 
function) versus length-scale. The scale r is normalized by the integral length L = 1.367. 
Different symbols are for different bins: o: hin\\ A: hin2\ ferns, O: hin^. The dashed lines 
are power-law fits yielding, in (a) for x = —0.3; in (b) for x = —0.49; and in (c) for 
IVpl, x = -0.55. 



4 Conclusions 



Using the JHU public turbulence database, we perform Lagrangian analysis of tempo- 
ral time correlations for the absolute values of velocity and pressure gradients. Con- 
sistent with earlier results for the square strain-rate and rotation rates |25ll26] . we find 
significantly longer decay times for the magnitudes as compared to the tensor or vec- 
tor elements, especially for strain-rate and pressure gradients. It is demonstrated that 
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Lagrangian dynamics of velocity gradient and pressure gradient (almost equivalent to 
acceleration since viscous effects are negligible) is determined mainly by scales provided 
by the locally averaged rate of dissipation, as predicated in the Kolmogorov Refined 
Similarity Hypothesis. We point out that a Lagrangian KRSH has also been shown 
to hold in the context of moments of two-time velocity increments [33]. That analysis 
is based on a fully Lagrangian rate of dissipation Et averaged over temporal domains 
of duration r along the particle trajectory. By its nature, er averages dissipation at 
various times. The present analysis is instead based on the more often used spatial 
average of dissipation at a single (initial-condition) time. Still, present results show 
that even after using local time-scales corrected for intermittency, there was remain- 
ing scatter observed in the correlation functions for absolute value variables. A better 
understanding of the origin of these deviations, as well as relating the relatively large 
intermittency corrections to various phenomenological models of turbulence, would be 
desirable developments. 
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